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ABSTRACT 

Observation of the H2O megamaser galaxy IC 2560 with the Chandra Ob- 
servatory reveals a complex spectrum composed of soft X-ray emission due to 
multi-temperature thermal plasma, and a hard continuum with strong emission 
lines. The continuum is most likely a Compton reflection (reprocessing) of pri- 
mary emission that is completely absorbed at least up to 7 keV. The lines can 
be identified with fluorescence from Si, S and Fe in the lowest ionization stages. 
The equivalent widths of the Si and S lines are broadly compatible with those 
anticipated for reprocessing by optically thick cold plasma of Solar abundances, 
while the large equivalent width of the Fe line requires some overabundance of 
iron. A contribution to the line from a transmitted component cannot be ruled 
out, but the limits on the strength of the Compton shoulder make it less likely. 
From the bolometric luminosity of the nuclear region, we infer that the source 
radiates at 1 - 10% of its Eddington luminosity, for an adopted central mass of 
3 X 10^ Mq. The overall spectrum is consistent with the hypotheses that the 
central engines powering the detected megamsers in accretion disks are obscured 
from direct view by the associated accretion disk material itself, and that there is 
a correlation between the occurrence of megamaser emission and Compton-thick 
absorption columns. For the 11 known galaxies with both column density mea- 
surements and maser emission believed to arise from accretion disks, eight AGN 
are Compton thick. 

Subject headings: masers — galaxies: active — galaxies: individual: IC 2560 — 
X-rays: galaxies 
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1. Introduction 

Position- velocity resolved interferometric mapping of water vapor maser emission {urest = 
22235.08 MHz) in active galactic nuclei (AGN) is a powerful tool for the study of gas struc- 
tures and dynamics in regions less than of order 1 pc from massive central engines. For 
the cases in which the megamaser emission can be identified with an accretion disk rather 
than with the wind- or jet-like outflow, such mapping may allow the determination central 
engine (and presumably black hole) mass, Eddington luminosity, and disk geometry (e.g., 
warping). The three most notable examples are: the geometrically thin accretion disk of 
NGC 4258 (Myioshi et al. 1995; Greenhill et al. 1995) which obeys a Keplerian rotation law 
to ^ 1%; the apparently massive accretion disk within the NGC 1068 AGN (Greenhill et 
al. 1997; Kumar 1999; Lodato & Bertin 2003); and the accretion disk within the nucleus of 
the Circinus Galaxy (Greenhill et al. 2003), and the associated wide angle wind. 

Water megamaser and nuclear activity are linked by the maser pump process which in 
the case of accretion disks, is plausibly driven by collisions in molecular gas heated by X-ray 
irradiation (Neufeld, Maloney, & Conger 1994). Nonetheless, detectable megamaser emission 
is relatively rare. This is probably a consequence of intrinsic luminosity and geometry. In 
general, maser emission is amplified most along the directions in which there is the most 
material with overlapping line-of-sight velocities, and as a result the emission is beamed 
narrowly along tangent planes for any individual accretion disk. Association with type-2 
AGN is a consequence. 

Average detection rates among type-2 AGN are on the order of 10% (e. g. Braatz et 
al. 2004), depending on galaxy distance and instrument sensitivity. Only ~ 60 megamaser 
systems are known (Kondratko et al. 2005b, c; Henkel et al. 2005 and references therein) 
and evidence from spectra and Very Long Baseline Interferometer (VLBI) maps associating 
the emission with accretion disks is available for about a third (e.g., Kondratko, Greenhill, 
& Moran 2005a). Nonetheless, the contribution to understanding of AGN by broadband 
studies of megamaser galaxies lies in the qualitatively greater level of detail that may be 
achieved, specifically because the structure and dynamics of the molecular gas can be so well 
constrained. 

Though the central engines associated with megamasers are likely to be obscured in the 
optical, UV, and soft X-ray bands - at least for the absorbing columns corresponding to less 
than a few Thomson optical depths - X-ray spectroscopy can provide a good estimate of the 
absorbing column, and thus the true luminosity of the central source. For other, more heavily 
obscured sources, the details of the profile of the Fe Ka line can often provide additional 
constraints (e.g., Levenson et al. 2002). One good candidate for detailed measurements is the 
Seyfert 2 galaxy IC 2560. It contains an H2O megamaser (Braatz, Wilson, & Henkel 1996), 
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with a peak flux density of up to ~ 0.4 Jy. This relatively nearby (26 Mpc; Vrec — 2876 km 
s~^) source has been partially resolved with the VLBI (Ishihara et al. 2001) and a central 
black hole mass of ~ 3 x 10® Mq inferred, which is relatively low but comparable to the mass 
of the black hole in our own Milky Way (Ghez et al. 2005; Schodel et al. 2003) and Circinus 
Galaxy (Greenhill et al. 2003). 

Prior to the observation described in Section 2. this object has been observed with 
Chandra in 2000 October, via a short (~ 10 ks) observation, reported in Iwasawa, Maloney, 
& Fabian (2002), and even before that, with Asca (Ishihara et al. 2001). The data reported 
here are based on a longer, 50 ks pointing obtained in 2004 February. 

2. Observations, data reduction, and spectral fitting 

A 50 ks Chandra observation of IC 2560 was conducted on 2004 February 16-17. The 
processed data were reduced using the XSELECT tool, with an independent analysis using 
ClAO. Data were collected with ACIS-S, with the source on chip S3. The intensity of 
unrejectcd background was constant and relatively low (less than 2 ct s^-*^ in the entire 
ACIS-S3 chip), showing no sign of any flares, so removal of any additional data segments 
was not necessary. The image consists mainly of the point source, but there is a hint of some 
extended component. The spectrum of the extended component is under investigation. To 
measure the spectrum of the nuclear source, we extracted the data from a region 2" (~ 250 
pc) in radius, and determined the background from an annulus with inner radius 30" and 
outer radius 90" centered on the source. This corresponds to roughly 3.8 kpc and 11.3 kpc 
at the source's distance, and while might also include some of the diffuse emission from the 
galaxy itself, in the quite compact source region chosen here, the contribution of such diffuse 
emission (or instrumental background) was negligible, less than 0.1% of the source count 
rate. Total net count rate in the 0.5 - 10.0 keV range was 0.025 ct s~^. The spectrum was 
subsequently rebinned to allow at least 20 counts in a new spectral bin. We note that there 
is no variability in the source flux during the Chandra observation. Furthermore, the source 
shows the same count rate - as well as flux, for the same assumed models - as reported in 
Iwasawa et al. (2002). We prepared the ACIS S3 redistribution matrices and the effective 
area files as appropriate for this observation. Since the data were taken relatively recently, 
the CTI correction has been applied in the course of the pipeline processing. 

The background-subtracted data are well-represented by a spectrum with three general 
components: a very hard continuum, a soft component with some soft X-ray emission lines, 
and an emission line complex between 6 and 7 keV, presumably due to Fe K-shell transitions. 
This is illustrated in Figure 1 where we assumed as a model a simple power law, absorbed by 



Table 1. Results of data modeling without photo- ionization 



model^ kTi^ (keV) kT2^ (keV) EW (keV)^ xVdof 



1 mk+po+ga 


'-'•'-''-'-0.09 


u.uiy_o.oi3 






-0.3 ±0.3 


2.7 ±0.5 


89.6/49 


2 mk+rf+ga 


0.581°;}° 


0.023lHf 






^•^^-0.39 


9 O-I-0.07 


90.2/49 


3 2mk+rf+ga 


0.65t°;J° 




n 14+0.05 

'-'•-'-^-0.056 


1^ 


2.8 ±0.2 


'J-'-'-o.e 


83.7/48 


4 2mk+rf+3ga^ 


0.651°;^^ 




14+0-04 

'-'•-'-^-0.063 




9 7-h0.2 
^- ' -0.3 


3.4 ±0.7 


71.1/44 


5 ab+2mk+rf+3ga+ds'* 
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^Those describe spectral forms incorporated into the model: mk - mekal, po - power law, 
ga - Gaussian line, rf - Compton reflection (pexrav), 3ga - three Gaussian lines, ds - Compton 
down-scattered line shoulder, ab - additional absorption at redshift 2; = 0.01 of host galaxy 

"^mekal model plasma temperature and metal abundances ' 

'^Power law photon spectral index, N{E) oc 

^EW is the equivalent width of the Fe 6.4 keV line 

^Parameter fixed; the abundances are not well constrained (0.2 < A < 1) 

^Three Gaussian lines are the Fe-K, Si-K, and S-K lines originating in neutral material, 
as described in the text 

^Parameter uncertainty pegged at lower limit 

''Host galaxy Nh = OMto, 2I X 10^^ cm ^; and upper limit to the downscattered/observed 
line intensity of ~ 0.15 



Table 2. Results of data modeling with photo-ionization models 
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1 xstar^ +rf+3ga 

2 xstar^ +po+rf+3ga 

3 xstar^ +mk+rf+3ga 

4 2xstar^ +rf+3ga 

5 2xstar^ +mk+rf+3ga 


1.6 
1.6 
3.0 
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1.4 
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2.4 ±0.5 
2.40 ±0.45 


152/49 
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^Spectral components: xstar - XSTAR photo-ionization model; other components as in 
Table 1 

'^og of ionization parameter in XSTAR model cn 

I 

'^Plasma temperature in mekal model 

'^Power law photon spectral index, N{E) oc 

^EW is the equivalent width of the Fe 6.4 keV line 

^Model parameterized by ionization parameter ^ and H column density, Nh- Fits show 
no dependence on Nh, uncertainties on log(^) are usually smaller than model grid spacing 
(Alog(^) = 0.2), and are not shown 
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the Galactic column of 6.5 x 10^° cm~^. The residuals clearly show the features above. This 
spectrum is similar to that reported by Iwasawa et al. (2002), but the longer exposure clearly 
reveals more details, with better resolved spectral features. A somewhat more complex and 
realistic model including a power law, collisionally ionized plasma (described as a mekal 
XSPEC model; see below) and Gaussian line, respectively, for the hard continuum, soft 
emission, and the Fe line, gives an energy power law index a = — 1.3 ± 0.3 (Table 1, Model 
1). This is extremely hard for the intrinsic X-ray spectrum of an AGN, but since IC 2560 is a 
megamaser source, the AGN might be obscured by a large column of absorbing material and 
the continuum measured here is probably not the primary spectrum. Specifically, the primary 
continuum might well be a power law, but we anticipate that it is heavily absorbed rather 
than particularly hard (as compared to unobscured Seyferts), and it may be accompanied 
by a Compton reflection component, which is likely to be present when the primary source 
photons are reprocessed by circumnuclcar material. If the absorption is particularly severe - 
~ 10^^ cm~^ or more for Solar abundances - the very hard spectrum measured above ~ 2.5 
keV could be just that due to the reflected component, although some additional flux might 
arise via transmission through the absorber. 

Below, we describe the hard continuum as pure Compton reflection. When we replace 
the power law component representing the primary continuum with a reflection component 
at a fixed 60° inclination (using the XSPEC model pexrav; see Magdziarz Sz Zdziarski 
1995), we obtain a — I.SILq^, which corresponds to a more reasonable intrinsic spectrum 
of the illuminating photons. The choice of a pure refiection model is also supported by the 
detection of the very strong, intrinsically narrow Fe K line (equivalent width ~ 2.5 keV, 
depending on the adopted continuum model, and cr < 40 eV) at an energy GAOtoi keV, 
which is consistent with neutral material and in agreement with Iwasawa et al. (2002). The 
soft component below 1 kcV is fairly well flt by a single temperature, collisionally ionized 
plasma with kT = 0.5 ±0.1 keV, but with unphysically low abundances A = 0.02 ± 0.01 
(Tible 1, Model 2). 

The strong iron line should also be accompanied by a noticeable Compton downscattered 
shoulder, forming a continuum centered at a rest frame energy of ~ 6.3 keV, with intrinsic 
width of ~ 0.14 keV (lUarionov et al. 1979). The strength of this continuum relative to the 
unscattered narrow fine core depends on details of the reflecting material (Matt 2002; see 
Fig. 2). This component has been seen in other refiection-dominated AGN such as Circinus 
Galaxy (Molendi, Bianchi, & Matt 2003), NGC 4945 (Done et al. 2003), and NGC 1068 
(Matt et al. 2004). We incorporated the shape of this downscattered continuum (lUarionov 
et al. 1979) as a local model in XSPEC, assuming a temperature of 1 eV for the reflecting 
material, and include it in the spectral flt in Table 1, Model 5. While this downscattered 
continuum is not significantly detected in our data, the 90% confidence upper limit on its 
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strength is ~ 15% of that of the narrow hne core. The low value indicates super-Solar iron 

abundance, with the observed flux arising from a reflector viewed at a high inclination (Matt 
2002), as may be anticipated for reflection from the far side of a flattened structure with 
fairly constant ratio of height to radius. 

In addition to the 6.4 keV hne, there is probably a marginally resolved hne at 6.7 keV 

as well, although it is not highly significant (and thus not included in any fits in Table 1): 
adding three parameters improves the best-fit only by 5. Close inspection of the unbinned 
data does not reveal whether the line is intrinsically broad or divisible into sub-components, 
but the latter is a more likely possibility. If real, the 6.7 kcV emission indicates very hot or 
highly ionized material, in contrast to the neutral material responsible for the 1.75 (Si) and 
2.32 keV (S) hues discussed below, and the Compton reflecting material, which has a lower 
inferred ionization parameter. Presence of 6.7 keV emission provides marginal evidence of a 
multi-temperature plasma in close proximity to the central engine. 

Rather than a very low abundance single-temperature plasma, the soft spectrum may 
be well flt by a model comprising nearly Solar abundance material, but also including sec- 
ond plasma component with kT ~ 0.1 keV. The new component contributes the additional 
continuum necessary to dilute model line strengths and match the observations with more 
typical abundances (Table 1 as Model 3). These two temperature components arc probably 
simply an approximation to the multi-temperature hot gas which is expected in starburst re- 
gions (e.g. Strickland & Stevens 2000). It should not be confused with the muIti-tcmpcrature 
plasma inferred from the Fc line emission, which originates on much smaller scales. We note 
here that even though a 2-temperature model is consistent with Solar abundances, it con- 
strains them poorly: spectral flts - where we require the same abundance of both mekal 
components - allow abundances ranging from 0.2 to 1. We note here that it is possible to 
fit the spectrum with a combination of thermal plasma and photo-ionized components, but 
at least some contribution from thermal plasma is required; we discuss this more extensively 
below. However, there is no combination of hot gas components which can account for the 
major residuals in the spectral fit at ~ 1.8 and 2.3 keV. This is because these features are at 
the energies expected from neutral Si and S, so cannot be produced by hot plasma emission. 
The obvious cool plasma component is the one which produces the refiected continuum. For 
Solar abundances this should produce an equivalent width of ~ 250 and 170 eV in neutral Si 
and S, respectively (Matt, Fabian, & Reynolds 1997) to accompany the Fe line of equivalent 
width 0.8-1.6 keV (depending on atomic data and assumed geometry: George & Fabian 
1991; Zycki & Czerny 1994, Matt, Brandt, & Fabian 1996). Adding in two additional Gaus- 
sian lines gives Si line emission with equivalent width of l^^tfo eV at 1.75 ± 0.03 keV and S 
of 200 ± 140 cV at 2.32 ± 0.12 keV, accompanying the neutral fiuorescent iron line of large 
equivalent width, 3.4 ±0.7 keV (Table 1, Model 4). Thus while the Si and S line energies 
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and strengths are consistent (though with fairly large error bars) with the fluorescent lines 
expected from reflection from X-ray illuminated neutral material with Solar abundances, 
the iron line is relatively strong, which indicates super-Solar abundances in this element. 

There is some evidence for additional cold photoelectric absorption, with the best flt 
yielding Nh ~ 50 - 70 x 10^° cm~^ in addition to the 6.5 x 10^° cm~^ attributable to our 
Galaxy. This inference depends on the description of the hot gas component: this is because 
the abundances and relative intensities in the mekal plasma are strongly correlated with the 
inferred level of photoelectric absorption in the spectral fitting. In any case, such columns 
are generally seen towards starbursts, as expected from the molecular gas in these regions 
(e.g., Pietsch et al. 2001). Note that this is in contrast to the primary continuum of the 
AGN which we assume to have very large absorption, > 10^^ cm~^. The model including 
such additional absorption is now given in Table 1 as Model 5. 

Although the observed soft X-ray emission is consistent with arising in a gas-rich, dusty 
starburst region, there is a potential alternative origin via scattering in partially ionized 
material. Such extended emission from photo-ionized gas is seen in the soft X-ray spectra 
of other AGN (e.g. Mkn 3; Sako et al. 2000). Here we investigate the origin of this emission 
by replacing the two mekal plasma components with the emission expected from photo- 
ionized gas (modeled using the XSTAR code: Bautista & Kallman 2001; gridl9c from the 
XSTAR web page). A simple model, with one XSTAR component replacing the two mekal 
components (from Model 4 in Table 1) gives a bad description of data, xl = 152/49 (Model 1 
in Table 2). The fit can be improved by adding a power law representing the primary emission 
scattered off a completely ionized plasma, but the fit is still unsatisfactory {xl = 106/48; 
Model 2 in Table 2), with strong residuals around 0.6 and 1 keV. However, previous fits 
indicate that the emitting plasma is multi-temperature, therefore we also try more complex 
models. Two photo-ionized XSTAR components (with different ionization parameters) give 
a better flt, xl — 96.4/47 (Model 4), but strong residuals still remain at low energies. These 
residuals cannot be explained by any pair of photo-ionized components, but an inclusion of 
an additional mekal component provides a satisfactory fit, giving xt = 58.1/45 (Model 5). 
For completeness, we have also investigated that a simpler hybrid model with one XSTAR 
and one mekal component (Model 3), but this gives an unsatisfactory fit, xl = 89.4/47. 
Thus, at least some contribution from mechanically-heated plasma is necessary to explain 
the soft X-ray emission, but even then, the photo-ionized plasma has to have a range of 
ionization stages. Presence of mechanically (coUisionally) heated plasma is in fact expected 
in IC 2560, as the galaxy contains a starburst region (Cid Fernandes et al. 2004). 

To quantify this, we investigated the relative fluxes of the mekal and photo-ionized 
plasmas. The two mekal components in Model 5 of Table 1 contribute comparably - at the 
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observed flux level of ~ 2 x 10"^^ erg cm~^ s~^. In the hybrid model, Model 5 of Table 2, 
the mekal component contributes comparably to that in low T component in the two-mekal 
model (~ 2 X 10"^^ erg cm~^ s^^), while the sum of the two XSTAR components contributes 
~ 5 X 10^^"^ erg cm^^ s^^. This leads us to the conclusion that the total contribution of 
photo-ionized emission to the soft X-ray flux of IC 2560 can be as little as 0% and much as 
~ 70% of the total, but (1) at least some coUisionally heated plasma is required, and (2) the 
photo-ionized component implies a wide range of ionization parameters. 

Hence we use the best fit starburst/refiection dominated AGN model (that given in 
Table 1 as Model 5) to plot the unfolded spectrum shown in Fig. 3. A remaining small 
residual can be seen at ~ 5 keV. Including a free Gaussian line to model this feature gives a 
better fit with xl = 46.9/39, formally a significant detection with line energy 5.1 ± 0.2 keV, 
width O.Sloi keV and very large EW of l-Sl^i^ keV. However, no plausible atomic features 
exist at this energy. At such high equivalent width it is most likely to be iron but the energy 
shift is very large, requiring the material to be flowing away from the observer at ~ 50, 000 
km s~^; a multiple Compton downscattered shoulder explanation seems even more contrived. 
However, this best fit also requires that the refiected emission is now due to an unusually 
steep power law illuminating the refiector, with energy spectral index a = 2.1. Assuming 
that the intrinsic AGN spectrum has a more usual a = 1 gives xl = 54.3/40, compared to 
57.7/43 without the additional fine, so it is not significant on an F-test. Thus it seems most 
likely that this feature may be an artifact in a complex, low signal-to-noise spectrum. 

With such a model, and including the emission lines as given below, the measured total 
2-10 keV fiux is 3.3 x 10"^^ erg cm^^ s^^. When the emission lines discussed below are 
accounted for separately, then measured reflection model flux (and this of course does not 
include the primary continuum) is 2.6 x 10~^^ erg cm~^ s~^. In the extreme case, if the 
reflector were to subtend a solid angle as large as 27r, it would have to be illuminated by an 
intrinsic, unabsorbed - but invisible directly to the observer - 2-10 keV flux of 4.5 x 10~^^ 
erg cm~^ s^^. A more likely geometry might be one in which much smaller fraction of the 
refiector is visible to us, in which case the unabsorbed fiux of the AGN is correspondingly 
higher. We can thus set a lower limit on the 2-10 keV luminosity of ~ 3 x 10**^ erg s~^: at 
least this much is needed to account for the refiected spectrum. 

3. Origin of the X-ray spectral features and the relationship of the megamaser 

emission to the X-ray absorption 

The spectral fits presented in the previous section indicate that the iron Ka line is some- 
what stronger than predicted, perhaps indicating super-Solar iron abundance. We quantify 



-10- 



this by a Monte Carlo simulation of reflection from a square cross-section torus-like struc- 
ture with Thomson optical depth tt = 3, with inner edge starting 50° from the "torus" axis 
(see Krolik, Madau, & Zycki 1994). This gives a predicted iron Ka line equivalent width 
of 1.5 keV (Anders-Grevesse abundances: Fe/H = 4.7 x 10^^), or 1.2 keV (Anders-Ebihara 
abundances: Fe/H = 3.3 x 10"^), for viewing angles close to the equatorial plane. To get a 
reflected iron line as large as that observed requires at least 2x Solar abundance of Fe. The 
large iron Ka line also should be accompanied by strong K/3 emission at 7.05 keV, whose 
strength is about 10% of that of the Ka line. While the signal-to-noise ratio at these high 
energies is not high, such a feature (with all parameters flxed or tied to Ka) reduces xt by 
3 to 57.8/42. 

An alternative explanation for a large iron line equivalent width is the contribution from 
a heavily absorbed transmitted spectrum as well as reflection (Ghisellini, Haardt, & Matt 
1994; Levenson et al. 2002). If the observed hue is indeed due to transmission, it would 
imply that the direct component, expected at 20-100 keV, should be quite strong, as in the 
case of NGC 4945 (Done, Madejski, & Smith 1996). However, with Solar abundances, this 
would require the model torus to subtend a rather large solid angle with a polar opening 
angle of < ±15° (Levenson et al. 2002). The Compton downscattered shoulder would also 
be stronger than that seen in reflected emission, up to 40% of the strength of the narrow 
line core (Matt 2002; see also Fig. 2), which is inconsistent with the 15% hmit on the 
downscattered continuum estimated here. 

From the spectral fits described in Tables 1 and 2, a picture emerges where the heavily 
absorbed primary continuum is reprocessed by nearly neutral material producing a Compton 
reflection component as well as fluorescence hues from Fe, S and Si. IC 2560 joins the rather 
smaU number of AGN (Mkn 3: Sako et al. 2000 and NGC 6552: Reynolds et al. 1994) 
in which the fluorescence lines of species with low atomic numbers (Si and S) are seen, 
and thus, reflected continuum in the 1-10 keV band dominates. The reflector is likely to 
be formed by the same optically thick material that obscures a direct view of the central 
engine. In IC 2560, H2O megamaser emission is also believed to arise from a close to edge-on 
structure (Ishihara et al. 2001), which by analogy to the well studied case of NGC 4258, 
is probably a parsec or sub-parsec radius annulus in an accretion disk. Geometrically, the 
simplest explanation of the reflector, absorber, and maser are that all three correspond to 
the same structure, the accretion disk. (See Greenhill et al. 2003 for discussion of similar 
circumstances in Circinus Galaxy.) We expect to see in reflection only a small "shver" of this 
X-ray illuminated material on the far side of the disc (a feature readily accommodated if the 
disk is warped), as has been observed in the megamaser galaxies NGC 4258 (Miyoshi et al. 
1995) and Circinus Galaxy (Greenhill et al. 2003). Alternatively, reflection could originate 
from low covering fraction, dense clouds above and below the obscuring maser structure that 
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blocks our view of the central engine. The former explanation has the advantage of being 
simpler: only one high column structure required, while the latter requires different regions 
of optically thick material. 

Regardless of the detailed location and geometry of the Compton reflector, the covering 
fraction of reflecting material is less than unity. This provides a lower limit to the intrinsic 
X-ray fiux, ~ 4 x 10^^^ erg cm~^ in the 2-10 keV band, corresponding to a minimum 
bolometric (0.01 - 100 keV) fiux of 6 x 10~^^ erg cm~^ and luminosity L = 5 x 10'^^ erg 
s~^. Adopting the central engine mass estimated by Ishihara et al. (2001), we can calculate 
L/LEdd > 0.01. If wc assume that all of the observed IR emission (L = 3 x 10^^ erg s~^) 
is reprocessed nuclear flux, then a very conservative upper limit is L/L^dd < 0.1. In reality 
this upper limit could be reduced to account for the IR flux that might arise in the circum- 
nuclear starburst region. With these constraints, we can infer that the accretion rate is more 
like that inferred for the common, broad-line Seyfert 1 galaxies (e.g. Czerny, Rozahska, & 
Kuraszkiewicz 2004) and the heavily obscured (Seyfert 2) megamaser galaxies NGC 3079 
(Kondratko et al. 2005b), NGC 4945 (GreenhiU et al. 1997), and NGC 1068 (Greenhfll & 
Gwinn 1997), than it is for the low L/L^dd core of the classic megamaser galaxy NGC 4258 
(e.g., Lasota et al. 1996). 

IC 2560 bears a startling resemblance to NGC 4945 (Done et al. 2003; Schurch, Roberts, 
& Warwick 2002) in showing a reflection dominated hard X-ray spectrum, together with 
starburst dominated soft X-ray emission. In IC 2560, there is evidence for a young (< 25 
MYr) stellar population (Cid Fcrnandcs et al. 2004), so it seems likely that there is hot gas 
from a starburst providing the majority of the soft X-ray flux, rather than photo-ionization. 
These characteristics differ from Mkn 3 (Sako et al. 2000) and NGC 1068 (Ogle et al. 2003), 
which have clear signatures of predominantly photo-ionized soft X-ray emission. In both 
cases, detection of polarized broad optical emission lines, indicates that there is a "mirror" 
directly reflecting UV and soft X-rays from the central source. Neither IC 2560 nor NGC 4945 
show broad polarized emission lines. It is worth noting here that both IC 2560 and NGC 4945 
are powered by relatively low-mass black holes and have modest nuclear luminosities. Low 
intensity in the photo-ionized component might be simply due to low luminosity of the 
central source. 

A plausible scenario for many megamaser galaxies (as well as other Seyfert 2s) is that 
we view the primary power law continuum of a ~ 1 through an absorber of appreciable 
optical thickness. This can provide an absorbing column due to the neutral gas that might 
be relatively small (10^^ cm~^), where the primary continuum shows a modest photoelectric 
absorption cutoff as in Mkn 3 or NGC 4258. Alternatively, the absorber can provide a 
column of several times 10^^ cm-^ or more, as seen in NGC 4945, NGC 1068, and NGC 3079 
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- and IC 2560 belongs to the latter group (see Table 3). For material that is optically thick 
to electron scattering, the primary continuum is depressed at all energies. For a column with 
Thomson optical depth of order unity {Nh = 3 — 10 x 10^^ cm~^), the attenuated primary 
spectrum can be seen above 10 keV (as in NGC 4945), but for higher columns it is completely 
suppressed (as in NGC 1068). The column in IC 2560 can only be constrained by the current 
Chandra data to be greater than or equal to that seen in NGC 4945, Nh > 3 x 10^"^ cm~^. 
Sensitive observations at higher X-ray energies are required in order to determine whether 
the primary emission can be detected above 10 keV or whether like NGC 1068, IC 2560 is 
completely obscured along the hne of sight. 

The possibility of a physical association between megamaser emission and obscuring 
material raises the question of whether there may be a correlation between the occurrence 
of visible megamaser emission and detection of large obscuring columns. Most importantly, 
maser emission is anisotropic, and for a thin, moderately warped disk, it would be beamed 
in a narrow solid angle about the tangent planes (contingent on favorably small gradients in 
line-of-sight velocity). As well, it is conceivable that maser amplification could be enhanced 
(i.e., longer gain paths or narrower beam angle) when large columns are present. Direct 
evidence that megamaser emission originates from disk material is available for about a third 
of the ~ 60 known masers, chiefly in the form of VLBl maps or spectra that display highly 
Doppler shifted line complexes symmetrically bracketing the systemic velocity {Vgys)- Of 
these, there are estimates of X-ray absorption columns for 11 galaxies. Most of these systems 
are Compton thick; eight exhibit absorption in excess of 10^^ cm~^ (Table 3). Overall, 
the distribution of column densities is substantially skewed with respect to the absorption 
distribution for typc-2 AGN in general (Bassani et al. 1999; Risaliti et al. 1999). However, 
three known Compton thin cases are notable: NGC 4258, NGC 4388, and NGC 4051. Each 
exhibits variation in column density (Fruscione et al. 2005; Elvis et al. 2004; McHardy 
et al. 1995). In the case of NGC 4258, the accretion disk has been shown to cross the 
line of sight to the central source and modest variability in Nh over months is believed to 
arise from dumpiness at radii comparable to that of the persistent maser emission, on the 
order of 10^-10^ Schwarzschild radii (Fruscione et al. 2005; Herrnstein et al. 2005). In 
contrast, NGC 4388 has displayed extreme variability in column density, on time scales as 
short as hours, and to values as low as 2 x 10^^ cm~^. The inferred radius of the absorber is 
< 100 Ksch (Elvis et al. 2004). The angular structure and persistence of the maser emission 
is not known, but comparison to that of the absorber would enable testing of whether both 
arise from the same (disk-like) structure. The third AGN, NGC 4051 is unusual in that it is a 
narrow line Seyfert 1 (NLSyl) yet hosts a maser source that has been suggested to arise in a 
central accretion disk (Hagiwara et al. 2003). Column density variation is evident (McHardy 
et al. 1995), but unhke NGC 4258 and NGC 4388, the absorber is ionized and unhkely to 
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lie physically close to the maser medium. It is unlikely that there is any edge-on structure 
comprising cold material that could support maser emission, and the characteristics of the 
maser spectrum is consistent with origin in a wind rather than a disk (see Table 3 note). 
This is notable given the importance of nuclear winds in the formation of NLSyl spectra. 
On the other hand, the absorber and maser could both trace different parts of a warped 
structure that is only moderately inclined at small radii and tangent to some lines of sight 
at large radii. A clumpy or frothy medium might also enable a line of sight to the central 
source and ionized absorber at the same time as long gain paths are achieved among clumps 
at larger radii. 



4. Summary and conclusions 

Chandra observation of the H2O megamaser galaxy IC 2560 reveals a complex X-ray 
spectrum. At low energies, it exhibits a soft spectral component that is unlikely to be 
characterized by a single-temperature, coUisionally-ionized plasma, as this would require 
anomalously low elemental abundances. Instead, the soft component probably arises from 
such plasma - presumably associated with the host galaxy - at a range of temperatures, 0.1 
keV < kT < 0.7 keV, with abundances close to Solar. It is also possible to fit the spectrum 
by a model that combines one coUisionally ionized plasma component and a complex photo- 
ionized scattering medium, with a range of ionization parameters. 

At higher energies, the spectrum is dominated by a very hard X-ray continuum, to- 
gether with several emission line features at energies consistent with neutral Fe, Si, and S 
Ka lines. The best interpretation of such a hard spectrum is that the primary continuum 
is entirely absorbed (at least in the Chandra bandpass), requiring Nh > 0.3 x 10^^ cm~^, 
and that we observe only the much harder Compton reflection of the hidden primary con- 
tinuum from optically thick, high column density material. Given that the H2O megamaser 
emission is believed to arise in a close to edge-on disk-like structure and that maser action 
requires a large column of material, the simplest interpretation is that all three high column, 
low-ionization structures (absorber, reflector and megamaser) arise from a single physical 
component, probably a pc- scale or smaller accretion disk. We hypothesize that this is the 
case more generally and anticipate a correlation of megamaser emission and high X-ray ob- 
scuration. In an initial sample of 11 galaxies, this hypothesis is borne out. Discovery of 
new maser sources, identification of emission origins for known sources, and measurement of 
X-ray spectra for host AGN will enable further testing. 

The Compton reflection from this material also produces strong fluorescent emission 
lines. The observed intensities of the Si and S lines are entirely consistent with reprocessing 
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Table 3. Disk masers in AGN with estimated column densities^^^ 



Galaxy 


ID (2) 




References 




VLBI 


Spectra 


(10^3 cm-2) 


Nh 


Maser 


NGC 4945 






45 ±4 


1 


13 


Circinus Galaxy 






431^ 


2 


14 


M 51 








3 


15 


NGC 1386 






> 22 


4 


16, 24 


NGC 1068 




%/ 


>ioo 


5 


17 


NGC 3079 






>ioo 


6 


18 


IC 2560 






>30 


7,8 


19 


NGC 3393 






44+25 


9 


20 


NGC 4258 






0.6-1.3 


10 


21, 22 


NGC 4388 






0.02 - 4.8 


11 


23 


NGC 4051 






(cold); 0.8 - 3.7 (ionized) 


12 


24 



(1)- Henkel et al. (2005) and Kondratko et al. (2005b, 2005c) list known H2O masers in AGN. 
Disk masers are those originating within accretion disks, at parsec or sub-parsec radii. 
^2^- Nature of evidence supporting disk maser identification. VLBI mapping provides direct 
evidence. Spectra exhibiting red and blucshifted emission (more or less) symmetrically bracketing 
Vsys provide indirect evidence. For inclusion here, offsets from systemic must be > ±100 kms~^. 

Notes on specific objects: 

NGC 4945- Disk-like structure detected, similar to NGC 3079. 

Circinus Calaxy- Column > 10^^ cm"^ is admitted in an alternate model. 

M 51- Red and blueshifted maser emission are intermittent separately. 

NGC 1386- VLBI detection of red side of the disk. Later detection of the blue side in spectra. 
NGC .^55^-Classified as a Seyfert 2 galaxy, but recently, exhibited large, rapid change in Nh- 
NGC 405 1- Narrow-line Seyfert 1; no neutral absorber; ionized absorber with significant variabil- 
ity. Spectrum comprises narrow lines distributed over ~ 280 kms~^, though not well centered on 
Vgys- Hagiwara et al. (2003) infer association with an accretion disk. However, broadly spread 
narrow lines are also consistent with emission from a wind, as in Circinus (Greenhill et al. 2003). 

References- (1) Madejski et al. (2000); (2) Matt et al. (1999); (3) Fukazawa et al. (2001); 
(4) Guainazzi et al. (2005); (5) Matt et al. (1997); (6) lyomoto et al. (2001); (7) Iwasawa et 
al. (2002); (8) this paper; (9) Guainazzi et al. (2005); (10) Fruscione et al. (2005); (11) Elvis 
et al. (2004); (12) McHardy et al. (1995); (13) Greenhill et al. (1997); (14) Greenhill et al. 
(2003); (15) Hagiwara et al. (2001); (16) Braatz et al. (1997); (17) Greenhill & Gwinn (1997); 
(18) Kondratko et al. (2005b); (19) Ishihara et al. (2001); (20) Kondratko et al. (2005c); (21) 
Miyoshi et al. (1995); (22) Greenhill et al. (1995); (23) Braatz et al. (2004); (24) Hagiwara et 
al. (2003) 
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by Solar abundance material, while the Fe K line, with equivalent width of ~ 2.5 keV, is a 
factor of 2 higher than expected. The line might arise in reflection, transmission, or both. 
However, the most likely explanation for the high equivalent width of the line is Compton 
reflection from neutral medium with super-Solar (~ 2x) Fe abundances. This is further 
supported by the estimated upper limit on the intensity of the Compton "shoulder" that 
should be associated with the Fe K line. 

The intensity of the reprocessed component provides a lower limit on the 2-10 keV flux 
of the unabsorbed primary continuum of 4x 10~^^ erg cm^^ s^^, corresponding to a lower limit 
on the 0.01 - 100 keV luminosity of ~ 5 x 10^^ erg s^^. Conversely the observed bolometric 
luminosity of ~ 3 x 10*^^ ers s~^ places an upper limit. Since the mass of the central source 
estimated from partial resolution of the megamaser angular structure is ~ 3 x 10^ Mq, we 
infer that the source is accreting at a moderate rate corresponding to 0.01 < L/Lem < 0.1, 
as is the case for other heavily absorbed megamaser galaxies (e.g., NGC 1068). Further clues 
to the structure of this source arc likely to be revealed by sensitive hard X-ray observations 
at 10 - 50 keV, with missions such as NASA's NuSTAR. 

This project was partially supported by Chandra grants no. G04-5125X and G04- 
5127X from NASA via Smithsonian Astrophysical Observatory, by the Polish KBN grants 
2P03D01225 and PBZ-KBN-054/P03/2001, and by the Department of Energy contract to 
SLAC no. DE-AC3-76SF00515. L. G. thanks SLAC as well as KIPAC (Stanford University) 
for the hospitality at the time when this research was conducted. 
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— Chandra count spectrum for IC 2560, fitted to a simple power law model, absorbed 
energies by the Galactic absorption of 6.5 x 10^*^ cm~^. The residuals to the fit (plotted 
contribution to x^) show strong residual features, modeled and discussed in Sec. 2. 
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Fig. 2. — Profiles of the Fe Ka line in different geometries, showing differences in the strength 
of the Compton shoulder. Solid (black) curve shows the profile from transmission through 
a vertically-extended structure (assumed to have a geometry of a torus), seen close to equa- 
torial plane, dominated by reflection from its opposite side. Red (dotted) is the profile for 
pure transmission through the "torus" in equatorial plane (opposite side obscured). Green 
(dashed) is the profile for usual Compton refiection seen at 30 degrees. Profiles have been 
normalized to equal fiuxes in the line core. 
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Fig. 3. — Unfolded spectrum of IC 2560, as collected by Chandra in 2004. The model 
spectrum includes the Compton reflection of the invisible continuum, the two thermal plasma 
components, and the emission lines as discussed in the text. 



